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^  The  flow  of  a  liquid  with  a  free  surface  over  a  weir  in  a  channel  is 
calculated  nxunerically  for  thin  weirs  in  channels  of  various  depths,  and  for 
broad  crested  weirs  in  channels  of  infinite  depth.  The  results  show  that  the 
upstre2un  velocity,  as  well  as  the  entire  flow,  are  determined  by  the  height  of 
the  free  surface  far  upstream  and  by  the  geometry  of  the  weir  and  channel,  in 
agreement  with  observation.  The  discharge  coefficient  is  computed  for  a  thin  - 
weir,  and  a  formula  for  it  is  given  which  applies  when  the  height  of  the  weir 
is  large  compared  to  the  height  of  the  upstream  free  surface  above  the  top  of 
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SIGNIFICANCE  AND  EXPLANATION 


The  flow  of  a  liquid  with  a  free  surface  over  a  weir  in  a  channel  is 
considered.  Observation  shows  that  the  upstream  velocity,  as  well  as  the 
entire  flow  are  determined  by  the  height  of  the  free  surface  far  upstream  and 
hy  the  geometry  of  the  weir  and  channel.  Therefore  the  discharge  of  fluid 
along  the  channel  can  be  found  without  measuring  the  upstream  velocity,  which 
is  one  of  the  reasons  for  using  weirs.  Op  to  the  present  this  rather 
surprising  fact  has  not  been  treated  hydrodynamically. 

To  analyze  it,  we  calculate  numerically  the  flows  over  thin  two 
dimensional  weirs  in  channels  of  different  depths  and  over  broad  crested  two 
filmensional  weirs  in  channels  of  infinite  depth.  The  numerical  results  show 
that  the  entire  flow  is  determined  by  the  height  of  the  free  surface  far 
upstream  in  agreement  with  observation.  The  discharge  coefficient  is  computed 
tor  a  thin  weir  and  a  formula  for  it  is  given  which  applies  when  the  height  of 
t-he  weir  is  large  compared  to  the  height  of  the  upstream  free  surface  above 
ihe  top  of  the  weir.  The  coefficients  in  this  formula  are  close  to  those 
fouiid  empiricallv. 


if'  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
nitnary  lies  with  MRC,  and  not  with  the  authors  of  this  report. 
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1.  Introduction 

The  flux  or  discharge  Q  of  fluid  along  a  channel  is  the  product  of  the  fluid’s  mean  velocity  U  and  cross- 
sectional  area  A{H),  which  depends  upon  the  free  surface  level  H: 

Q  =  UA(H).  (1.1) 

Both  U  and  H  must  be  known  to  determine  Q  from  (1.1).  However  observation  shows  that  in  a  channel 
partially  obstructed  by  a  weir,  Q  is  determined  by  H  alone.  Then  Q  can  be  found  without  measuring  U, 
which  is  one  of  the  reasons  for  using  weirs  [Ackers  et  al.  (1978)]. 

This  observation  implies  that  although  U  and  H  are  independent  in  an  unobstructed  channel,  they  are 
not  independent  in  a  channel  containing  a  weir.  Instead  U  =  U{H)  is  a  function  of  H  and  of  the  size  and 
shape  of  the  weir.  Up  to  the  present,  this  rather  surprising  fact  has  not  been  treated  hydrodynamically. 

To  analyze  it,  we  shall  calculate  numerically  the  flows  over  thin  two  dimensional  weirs  in  channels  of 
different  depths,  and  over  broad  crested  two  dimensional  weirs  in  channels  of  infinite  depth.  This  requires 
solving  free  surface  problems  in  the  presence  of  gravity.  For  this  purpose  we  assume  that  the  fluid  is 
incompressible  and  inviscid,  and  that  each  flow  is  irrotational,  steady  and  two  dimensional.  We  also  assume 
that  there  are  no  surface  waves. 

For  a  thin  weir  at  each  surface  height  H,  we  find  that  there  is  a  solution  for  only  one  value  of  U,  which 
depends  upon  H  and  upon  the  height  W  of  the  weir.  This  confirms  the  empirical  result  mentioned  above. 
We  determine  the  flux  Q,  the  discharge  coefficient  C,  &nd  the  shape  of  the  free  surface  for  this  solution.  We 
also  obtain  a  formula  for  C  when  H/W  is  small,  where  H  is  measured  from  the  top  of  the  weir.  This  formula 
is  of  the  form  used  in  practice,  and  the  computed  coefficients  in  it  ru’e  close  to  the  empirical  ones  [Ackers  et 
al.  (1978)]. 

For  a  broad  crested  weir  of  length  L  and  height  W  =  oo,  we  find  that  there  is  a  unique  flow  for  given 
values  of  H  and  L.  We  determine  Q,  C  and  the  free  surface  for  two  values  of  H/L. 
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In  Section  2  we  define  C  and  introduce  the  form  of  the  equation  for  it.  Then  in  Section  3  we  consider  a 
thin  weir  in  a  channel  of  infinite  depth.  After  formulating  the  flow  problem  mathematically,  we  present  our 
numerical  procedure  for  solving  it  and  display  some  of  the  numerical  results.  In  Section  4  we  treat  a  thin 
weir  in  a  chemnel  of  finite  depth.  In  Section  5  we  treat  broad  crested  weirs  in  a  channel  of  infinite  depth. 
Finally  in  Section  6  we  discuss  a  related  free  surface  flow  that  can  be  analyzed  explicitly. 

2.  Discharge  coefficient 

Figure  la  shows  a  side  view  of  a  channel  with  a  thin  weir,  and  Figure  lb  shows  a  cross-section  of  the  channel 
at  the  weir.  The  height  H  of  the  free  surface  far  upstream  and  the  depth  W  of  the  bottom  are  measured 
from  some  point  on  the  upper  edge  of  the  weir,  as  is  shown  in  Figure  la.  The  cross-sectional  area  A(H)  of 
the  channel,  and  the  area  a(H)  of  the  region  above  the  weir  up  to  the  level  H,  are  shown  in  Figure  lb. 

In  terms  of  these  quantities  and  the  acceleration  of  gravity  g,  we  define  the  dimensionless  discharge 
coefficient  C  by 

Q  =  CigH)^a{H).  (2.1) 

The  fact  that  the  flow  is  determined  by  H  means  that  C  is  a  function  only  of  H  and  of  the  geometry  of 
the  weir  and  the  channel.  Since  C  is  dimensionless,  we  shall  write  it  as  a  function  of  the  dimensionless  ratio 
H/W,  C  =  C{H/W).  It  also  depends  upon  other  dimensionless  geometrical  parameters. 

We  assume  that  C  is  regular  for  H/W  small,  so  that  it  has  the  form 

C{H/W)  =  C(0)  -1-  C\G)HIW  -I-  0[{HIWf].  (2.2) 

For  flows  over  two  dimensional  thin  weirs  this  assumption  is  confirmed  by  our  numerical  results  ,  which  are 
presented  in  Sections  3  and  4.  The  results  also  determine  the  coefficients  C(0)  and  C'(0)  for  those  particular 
weirs. 

Let  us  specialize  these  formulas  to  a  weir  with  a  rectangular  opening  of  width  b,  for  which  a{H)  =  bH. 
Then  (2.1)  and  (2.2)  yield 

Q  =  g^Hh{C{0)  +  C'{0)H/W  +  O[iH/wy]}.  (2.3) 

For  a  weir  with  a  V-notch  of  angle  6,  a{H)  =  tan  5/2  and  we  get  instead 

Q  =  tan 5/2{Cv(0)  +  C'v{0)H/W  +  0[{H/Wf]}  (2.4) 

Results  of  the  forms  (2.3)  and  (2.4)  are  used  in  practice  [Ackers  et  al.  (1978)].  They  are  derived  by 
dimensional  analysis,  hydraulic  approximations  and  fitting  to  experimental  data.  Corrections  for  viscous 
and  surface  tension  effects  are  often  included. 

By  combining  (1.1)  and  (2.1),  we  can  determine  the  velocity  U.  From  it  we  find  that  the  Froude  number 
F  =  U[g(H  +  lV’)]"i  of  the  flow  in  the  channel,  based  upon  the  total  depth,  is 
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Figure  la:  Side  view  of  a  channel  with  a  thin  weir. 


Figure  lb: 

Cross-section  of  a  rectangular  channel  of  width  B 
at  the  weir,  which  has  rectangular  opening  of  width 
b,  at  height  W  above  the  bottom.  The  water  depth 
far  upstram  Is  W  +  H. 


F  =  C{H/W){1  +  W/H)-ia{H)/A(H) 

=  ^  +  0[{H/W)^}. 

For  a  weir  with  a  rectangular  opening  of  width  6  in  a  rectangular  channel  of  width  B,  we  have  a(ff)  =  bH, 
A{H)  =  B{H  +  W)  and  (2i5)  becomes 

^  (2.6) 

For  a  thin  weir,  the  only  other  geometrical  ratios  upon  which  C(0)  and  C'(0)  depend  in  this  case  are  b/B 
rmd  the  lateral  position  of  the  opening.  This  result  shows  that  F  is  small  of  order  (/f/W)i  when  H/W 
is  small,  so  then  the  flow  is  extremely  subcritical.  When  b  =  B  the  flow  is  two  dimensional  and  then  the 
coefficients  are  constants.  We  shall  determine  them  in  Sections  3  and  4. 

3.  Thin  weir  in  infinitely  deep  channel 

Let  us  consider  a  thin  weir,  such  as  that  shown  in  Figure  la,  in  a  channel  of  infinite  depth,  so  that  W  =  oo. 
Let  the  channel  be  rectangular  with  width  B  =  b,  eo  that  the  opening  above  the  weir  extends  completely 
across  the  channel.  Then  the  flow  is  two  dimensional.  We  seek  a  flow  which  becomes  a  falling  jet  with  two 
free  surfaces  after  it  crosses  the  weir.  Thus  the  flow  appears  as  in  Figure  la  with  the  bottom  removed  to 
infinity,  and  with  the  jet  falling  to  infinity. 

We  introduce  cartesian  coordinates  with  the  x-axis  directed  vertically  downwards  through  the  separation 
point  5,  and  with  the  asymptote  to  the  upper  freee  surfac^as  the  ^axis  (see  Figure  la).  Gravity  acts  in 
the  x-direction.  As  y  -*  -oo,  the  velocity  approaches  zero.  As  y  — ►  +oo  the  flow  approaches  the  thin  jet 
solutions  of  Keller  and  Weitz  (1957)  and  Keller  and  Geer  (1973). 

Let  the  potential  function  be  ip  and  the  stream  function  be  Without  loss  of  generality  we  choose 
tl>  =  0  on  the  lower  free  surface  and  ^  =  0  at  the  separation  point  5.  Let  Q  be  the  value  of  ip  on  the  upper 
free  surface.  On  the  two  free  surfaces  the  Bernoulli  equation  yields 

i{Vp)^-gx  =  0  (3.1) 

We  introduce  dimensionless  variables  by  taking  (Q*/y)^  as  the  unit  length  and  (Qs)i  as  the  unit 
velocity.  In  these  new  variables  (3.1)  becomes 

(V^)*-2x  =  0  on  V*  =  1  and  on  V*  =  0,  ip>0  (3.2) 

The  plane  of  the  dimensionless  potential  f  =  p  +  iip  is  shown  in  Figure  2. 

Let  the  complex  velocity  be  (^  =  u  —  it>.  Here  u  and  w  are  the  x  and  y  components  of  the  vector  velocity. 
As  /  -»  -00,  we  require  that  there  be  no  waves,  so  the  velocity  C  vanishes  like  .  As  f  —  +oo,  the  velocity 
C  increases  like  [Keller  and  Weitz  (1957)].  Thus  we  have 


Figure  2: 

The  complex  potential  plane.  For  a  thin  weir  In  an  Infinitely  deep 
channel,  the  point  E  coincides  with  I  at  -®.  For  a  channel  of 
finite  depth,  E  Is  on  the  axis  between  I  and  S,  as  It  Is  also 
for  a  broad  crested  weir  In  a  channel  of  Infinite  depth. 


^  ~  e7 

as  f  —*  — oo 

(3.3) 

as  /  — ♦  -boo 

(3.4) 

The  problem  is  to  find  ^  as  an  analytic  function  of  f  =  irl>  in  the  strip  0  <  V*  <  1.  satisfying  (3.2)-(3.4) 
and  the  kinematic  condition 


t/  =  0  on  V’  =  0,  (f  <  0 


(3.5) 


We  define  the  new  variable  t  by  the  relation 


/  = 


2(t2  + 1) 


(3.6) 


The  transformation  (3.6)  maps  the  flow  domain  into  the  interior  of  the  unit  circle  in  the  (-plane  so  that  the 
vertical  wall  goes  onto  the  real  diameter  and  the  free  surface  goes  onto  that  portion  of  the  circumference 
lying  in  the  upper  (-plane  (see  Figure  3). 

Following  deBoor  (1961)  we  define  the  function  n(()  by  the  relation 


<  =  -((-H)(-lnc(l-K*)]ie«(‘) 


(3.7) 


Here  c  is  a  real  constant  between  0  and  j.  We  shall  choose  c  =  0.2.  It  can  be  checked  eeisily  that  the  expression 
(3.7)  satisfies  the  conditions  (3.3)  and  (3.4).  The  function  Q{t)  is  analytic  for  |((  <  1  and  continuous  for 
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Figure  3; 

The  complex  t-plane.  For  a  thin  weir  in  an  infinitely  deep 
channel,  E  coincides  with  I  at  t  -  -1.  For  a  channel 
of  finite  depth  E  is  on  the  diameter  between  I  and  S. 
The  same  is  true  for  a  broad  crested  weir  in  a  channel  of 
infinite  depth. 


|<|  <  1.  The  kinematic  condition  (3.5)  implies  that  the  expansion  of  fl(t)  in  powers  of  t  has  real  coefficients 
With  this  expansion  inserted,  (3.7)  becomes 

c  =  -(t  +  1)[-  lnc(l  +  t»)]i  exp(f;  U„n  (3.8) 

n=0 

The  function  (3.8)  satisfies  (3.3)-(3.5).  The  unknown  real  coefRents  U„  have  to  be  determined  to  make  (3.8) 
satisfy  the  Bernoulli  condition  (3.2). 

We  use  the  notation  t  =  Iflc*"  so  that  points  on  the  free  surfaces  are  given  by  t  —  e"’ ,  0  <  ff  <  w.  Using 
(3.6)  and  the  identity 

we  obtain  after  some  algebra 

dx  _  I  si'n<r  «(<t) 

dff  2t  cos<rcos*  I  u(o-)^  +  t5((T)2  (3.10) 

dy  _  \  sin(T  6(<t) 

dff  2x  cos<TC08*  j  ti(<7-)2  +  t;((r)2  (3.11) 

Here  C(<^)  =  «(<t)  -  iv(ff)  denotes  the  value  of  <  at  a  point  on  a  free  surface. 


■■  -■  V.  ■' 
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The  value  of  x(<r)  on  the  free  surface  is  obtained  by  int<;grating  (3.10): 


x(a)  =  X,  + 


1  sinr _ v{t)  ^ 

2*  cos  T  cos*  5  «(»■)*  +  »(*)* 


0<<r<l 


(3.12) 


sinr  v{r)  j  x 

cos  T  cos*  ^  ti(T)*  +  €(z)*  2  ^  ^ 


(3.13) 


Here  x,  is  the  value  of  x  at  the  separation  point  5. 

We  solve  for  the  Un  numerically  by  truncating  the  infinite  series  in  (3.8)  after  N  terms.  For  convenience 
we  choose  N  to  be  even.  To  get  equations  for  the  coefficients  Un  we  use  collocation.  Thus  we  introduce  the 
N  mesh  points 

=  ^  +  ^=1 . N 


Using  (3.8),  (3.12)  and  (3.13)  we  obtain  and  x(<t/)  in  terms  of  the  coefficients  Un  nnd  x,.  Substituting 
these  expressions  into  (3.2)  at  the  point  <r/,  we  obtain  N  nonlinear  algebraic  equations  for  the  lV+1  unknowns 
UiiU], . Us-  Another  equation  is  obtained  by  integrating  along  the  equipotential  ^  =  0  from  V*  =  0 
to  V*  =  1  and  equating  the  value  of  x  at  V*  =  1  to  the  corresponding  value  of  x((t)  obtained  from  (3.12). 

We  solve  this  system  by  Newton’s  method.  Once  it  is  solved,  we  obtain  the  shape  of  the  free  surfaces 
in  parametric  form  by  integrating  numerically  the  identities  (3.10)  and  (3.11).  A  typical  profile  obtained  in 
this  way  is  shown  in  Figure  4. 

To  determine  the  coefficient  C(0),  which  occurs  in  (2.3),  we  calculate  x,,  the  x  coordinate  of  the 
separation  point  5.  It  follows  from  (2.3)  and  our  choice  of  dimensionless  variables  that  C(0)  is  related  to  x, 
by 

C(0)  =  x;K  (3.14) 


By  using  (3.14)  and  the  calculated  values  of  x,,  we  have  computed  C(0)  for  different  numbers  of  mesh  points 
N.  Figure  5  is  a  graph  of  the  values  of  C(0)  versus  1/N.  It  indicates  that  C(0)  varies  linearly  with  l/N  for 
N  large.  The  broken  line  in  Figure  5  represents  a  linear  extrapolation  to  =  oo,  which  gives 


C(0)  -  0.583 


(3.15) 


The  value  of  C'(0)  will  be  obtained  in  the  next  section. 


4.  Thin  weir  in  channel  of  finite  depth 

We  now  generalize  the  procedure  of  Section  3  to  include  the  effect  of  finite  depth,  so  that  W  oo  (see  Figure 
la).  The  dimensionless  potential  plane  is  shown  in  Figure  2.  By  using  the  transformation  (3.6),  we  map  the 
flow  domain  into  the  interior  of  the  unit  circle  in  the  (-plane,  so  that  the  horizontal  bottom  and  the  vertical 
wall  go  repectively  onto  the  portions  IE  and  ES  of  the  real  diameter.  The  free  surface  goes  onto  the  portion 
of  the  circumference  lying  in  the  upper  (-plane  (See  Figure  3). 


y 


Figure  4: 

«  «fnfile  for  the  flow  past  a  thin 
an'r„frnu:r,  U  c.«™el;  The  vertical  scale 
is  the  same  as  the  horizontal  scale. 


.  r  —  -en  the  complex  velocity 

Section  3  with  the  relation  (3.7)  replaced  by 


A.  In  S«tion  3  «.  cho«  c  -  0.2.  «,vi,e.  that  Ih.re  be  no  wave,  on 

..riiriUroircb-ollbel— eevabab.^ 

u.  Ih.  ..me  <(-l)  ot<  at  <  =  -1  by  the  .el.l.on 
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Figure  5; 

Numerical  values  of  C(0)  versus  1/N  obtained  by  the  scheme  of 
Section  3.  The  broken  line  corresponds  to  a  linear  extrapolation 
to  N  •  *. 


=  (4.2) 

The  procedure  of  Section  3  yields  N  +  1  equations  for  the  N  +  2  unknown  x,,eo,Ui, . . .  ,Un-  Another 
equation  is  obtained  by  using  (4.2)  where  the  Froude  number  F  is  specified.  A  typical  profile,  for  F  =  0.1,  is 
shown  in  Figure  6.  The  corresponding  value  of  H/W  is  found  to  be  0.42.  As  F  — ►  oo  we  find  that  eo  — ►  —1 
and  the  solutions  approach  the  solution  presented  in  Section  3. 

In  order  to  check  the  formula  (2.6)  and  to  determine  the  coefficients  C(0)  and  C'(0),  we  plotted 
versus  H/W  in  Figure  7  for  small  values  of  H/W.  These  values  were  obtained  by  using 
the  extrapolation  procedure  of  Figure  5.  The  curve  in  Figure  7  is  very  close  to  a  straight  line  of  slope  —0.80 
which  intersects  the  vertical  axis  at  0.583.  Therefore 

C(0)~  0.583  (4.3) 

C'(0)  -  |C(0)  ~  -0.80  (4.4) 


The  value  of  C(0)  predicted  by  (4.3)  agrees  with  the  value  (3.15)  obtained  in  Section  3.  From  (4.3)  and  (4.4) 
we  find 


Figure  6: 

Computed  free  surface  profile  for  the  flow  past  a  thin  weir  In 
a  channel  of  finite  depth  for  F  »  0.1.  The  vertical  scale  Is 
the  same  as  the  horizontal  scale. 


C'(0)  -  0.07  (4.5) 

Ackers  et  al.  (1978,  p.  57)  present  experimentsJ  values  of  C  obtained  by  various  investigators  for 

H/W  =  0  and  for  H/W  =  1.  As  (2.2)  shows,  these  are  just  the  values  of  C(0)  and  C(0)  +  C'(0)  H - .  From 

the  values  listed  in  their  table  (3.1)  we  find  that  the  experimental  values  of  C(0)  and  C'(0)  are  in  the  ranges 

0.564  <  C(0)  <  0.591  (4.6) 


0.066  <  C'(0)  <  0.085 


Our  numerical  results  (4.3)  and  (4.5)  are  evidently  within  the  experimental  ranges. 


(4.7) 


F/(H/W) 


Figure  7;  Values  of  F(H/W) 


versus  H/W. 


5.  Broad  crested  weir  in  water  of  infinite  depth 

We  shall  now  extend  the  procedures  of  Sections  3  and  4  to  obtain  the  flow  past  a  broad  crested  weir  in 
water  of  infinite  depth  (See  Figure  8).  We  denote  by  L  the  dimensionless  length  of  the  weir  and  we  write 
the  flux  Q  in  the  form  (2.1)  with  the  discharge  coefficient  Cb{H/L).  The  dimensionless  potential  plane  and 
the  complex  t-plane  are  the  same  as  in  Section  4  (See  Figures  2  and  3). 

As  t  — ►  -eo,  the  complex  velocity  ^  grows  like  (t  +  eo)"'/*.  As  <  — ►  —  1,  C  vanishes  like  t  +  1.  Therefore 
we  replace  relation  (4.1)  by 

a  =  i[-  In  c(l  +  t»))‘/’(t  +  1)(<  +  eo)-‘'*e"(‘)  (5.1) 


We  then  follow  step  by  step  the  procedure  of  Section  4  with  (4.2)  replaced  by  an  equation  which  expresses 
the  fact  that  the  distance  between  E  and  5  is  equal  to  L.  This  equation  is  found  by  integrating  numerically 
the  identity  (3.9)  along  the  streamline  V*  =  0  between  E  and  5.  Then  Cb  is  found  from  (3.14). 

A  typical  profile  for  H/L  =  0.88  is  shown  in  Figure  8.  A  similar  profile  was  found  for  H/L  =  1.23.  The 
corresponding  values  of  the  discharge  coefficient  are 


Figure  8: 

Computed  free  surface  profile  for  the  flow  past  a  broad 
crested  weir  for  L  »  1.6.  The  vertical  scale  Is  the  same 
as  the  horizontal  scale. 


Cb(0.88)  =  0.583,  (5.2) 

Cb(1.23)  =  0.617.  (5.3) 

These  results  indicate  that  for  broad  crested  weirs,  the  discharge  Q  depends  only  on  H  and  on  the  geometry 
of  the  weir  and  of  the  channel. 

6.  Discussion 


The  weir  flows  which  we  have  found  in  Sections  3-5  are  subcritical  free  surface  flows  without  waves.  The 
condition  that  there  be  no  waves,  which  we  imposed  in  our  formulation,  was  essential  in  determining  a 


particular  solution.  Therefore  it  is  of  interest  to  consider  a  case  which  can  be  treated  exactly,  and  which 
leads  to  results  of  the  form  (2.2)  and  (2.6). 

We  consider  the  flow  in  Figure  9,  which  is  bounded  by  one  free  streamline  and  two  horizontal  rigid  walls. 
We  assume  that  there  are  no  waves  and  we  integrate  the  x-component  of  the  Euler  equation  of  motion  over 
the  flow  domain.  Upon  using  the  divergence  theorem  we  obtain 

^[«ntix  +  (/>"‘p+ffy)n,]ds  =  0.  (6.1) 

This  line  integral  around  the  boundary  of  the  flow  domain  can  be  evaluated  as  follows.  First  the  normal 
velocity  u„  vanishes  on  the  walls  and  on  the  free  streamline,  while  u„  =  Ur  =  Ui  or  U2  on  the  vertical  lines 
at  X  =  ±00.  Next  we  see  that  the  x-component  of  the  unit  normal  vanishes  on  the  walls  while  yn^ds  =  ydy 
can  be  integrated  along  the  free  streamline  where  p  =  0.  Finally  p~^p  +  gy  =  y(/fj  -  Hi)  at  x  =  +00  and 
=  g{H2  —  Hi)  +  (U|  —  Ui)/2  at  X  =  -00  from  the  Bernoulli  equation.  By  using  all  these  facts  in  (6.1)  we 
And 

-U^Hi  -  [(Ui  -  Ui)/2  +  g(Hi  -  /f j)]/f,  -  s(/fj  -  Hi)^/2  +  +  g(H2  -  Hi)H2  =  0.  (6.2) 


Figure  9: 

Sketch  of  a  two-dimensional  flow  emerging  from  beneath  a  flat  plate 
above  a  horizontal  bottom. 


Now  we  combine  (6.2)  with  the  mass  conservation  equation  UiHi  =  U2H2  to  get 

iH2-Hi)H9Hi-Ui)  =  0.  (6.3) 

Thus  either  H2  =  Hi  or  U2  =  (yffi)i.  In  the  second  case  we  see  that  the  flow  velocity  is  determined  by  the 
flow  geometry.  Then  the  Froude  number  F  =  U2(gH2)~^  and  the  flux  Q  =  U2H2  are  also  determined  by  the 
flow  geometry.  To  exhibit  the  similarity  of  the  results  to  those  for  a  weir  we  set  Hi  =  W  and  //j  =  W  +  H. 
Then  we  find 


F  =  (l  +  i//H^)-i  =  l-^  +  0(^)*,  (6.4) 

Q  =  =  (6.5) 

When  H  >  0  this  flow  is  subcritical  and  without  waves,  just  like  the  weir  flows  we  have  calculated.  However 
for  H  <  0,  the  present  flow  is  supercritical. 

Vanden-Broeck  (1980)  has  included  waves  in  this  problem  in  the  particular  case  W  =  oo,  and  calculated 
their  amplitude  by  using  this  method.  Finally  we  note  that  uniqueness  theorems  for  weir  flows  are  lacking. 
For  flows  under  sluice  gates  some  uniqueness  results  have  been  obtained  by  Budden  and  Norbury  (1982). 
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